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ABSTRACT 

We illuminate dynamical properties of Kuiper Belt Objects (KBOs) in the 
3:2 ("Plutino") and 2:1 ("Twotino") Neptunian resonances within the model of 



resonant capture and migration. We analyze a series of numerical integrations, 
each involving the 4 migratory giant planets and 400 test particles distributed 
throughout trans-Neptunian space, to measure the efficiencies of capture into 
O ! the 3:2 and 2:1 resonances, the efficiencies of capture into Kozai-type secular 

resonances, and the libration centers and amplitudes of resonant particles, all 
as functions of the migration speed. We synthesize instantaneous snapshots of 
the spatial distribution of ~10 4 resonant KBOs, from which we derive the lon- 
gitudinal variation of the sky density of each resonant family. Twotinos cluster 
±75° away from Neptune's longitude, while Plutinos cluster ±90° away. Such 
longitudinal clustering persists even for surveys that are not volume-limited in 
their ability to detect resonant KBOs. Remarkably, between —90° and —60° of 
Neptune's longitude, we find the ratio of sky densities of Twotinos to Plutinos to 
be nearly unity despite the greater average distance of Twotinos, assuming the 
two resonant populations are equal in number and share the same size, albedo, 
and inclination distributions. We couple our findings to observations to crudely 
estimate that the intrinsic Twotino population is within a factor of ~3 of the 
Plutino population. Most strikingly, the migration model predicts a possible 
asymmetry in the spatial distribution of Twotinos: more Twotinos may lie at 
longitudes behind that of Neptune than ahead of it. The magnitude of the asym- 
metry amplifies dramatically with faster rates of migration and can be as large 
as ~300%. A differential measurement of the sky density of 2:1 resonant objects 
behind of and in front of Neptune's longitude would powerfully constrain the 
migration history of that planet. 
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Subject headings: Kuiper Belt — comets: general — minor planets, asteroids - 
celestial mechanics 



1. INTRODUCTION 

The substantial eccentricity, e p , of Pluto's orbit is well explained by Malhotra's (1995) 
theory of resonant capture by Neptune. In this scenario, Neptune migrated radially out- 
wards from the sun by scattering planetesimals towards Jupiter, captured Pluto into its 
3:2 mean-motion resonance, and amplified e p upon continuing its migration. The resonant 
amplification of e p can be understood either mechanistically using Gauss's equations (see, 
e.g., Peale 1986), or in terms of the preservation of an adiabatic invariant (see, e.g., Yu & 
Tremaine 2001). The discovery of dozens of Kuiper Belt Objects (KBOs) that share the 
3:2 Neptunian resonance with Pluto and that also exhibit large orbital eccentricities (Jewitt 
& Luu 2000) apparently vindicates this proposal that Neptune plowed its way outwards 
through a field of planetesimals early in the history of the solar system (Fernandez & Ip 
1984). 

In the particular numerical simulation presented by Malhotra (1995), the 2:1 Neptunian 
resonance is predicted to be about as equally populated as the 3:2 resonance. For brevity, we 
will refer to KBOs in the latter resonance as "Plutinos" and KBOs in the former resonance 
as "Twotinos." As of July 4, 2002, the Minor Planet Center (MPC) database contains 
~41 Plutino candidates (objects observed at multiple oppositions having fitted semi-major 
axes, a, within 0.3 AU of the exact 3:2 resonance location at = 39.5 AU), and ~5 
Twotino candidates (objects whose a's lie within 0.3 AU of the exact 2:1 resonance location at 
^2:1 = 47.8 AU). 1 The Twotino candidates possess substantial orbital eccentricities, e ~ 0.2- 
0.4, in accord with the predictions of resonant capture and migration. Membership in a 
resonance is confirmed by verifying that the appropriate resonant argument librates rather 
than circulates (see §2); an orbit classification scheme based on this more rigorous criterion is 
currently being developed by the Deep Ecliptic Survey (DES) Team (see Millis et al. 2002). 
For the present paper, we will consider the observed Twotinos to be outnumbered by the 
observed Plutinos by a factor of F 5 S ~ 8. Part of this bias must simply reflect the fact that 
a 2-.i > 03:2; & U other factors being equal, more distant objects are fainter and more difficult 
to detect. But part of this bias may also reflect selection effects that depend on the longitude 



1 http://cfa- www. harvard.edu/iau/lists/TNOs. html 
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and latitude of observation. A mean-motion resonant object will be preferentially found at 
certain locations with respect to Neptune — what we will call "sweet spots" on the sky. The 
sweet spots for Twotinos are not necessarily those of Plutinos. 

Ida et al. (2000) point out that Neptune's ability to resonantly capture objects varies 
with migration timescale. A migration timescale that is 20 times shorter than that consid- 
ered by Malhotra (1995) is found to severely reduce the probability of capture into the 2:1 
resonance. For the 3:2 resonance, the capture probability is affected less dramatically by 
reductions in the migration timescale. The relative robustness of the 3:2 resonance com- 
pared to the 2:1 resonance is explored analytically by Friedland (2001), who underscores the 
importance of the indirect potential for the latter resonance. 

This paper quantifies, within the confines of the model of resonant migration, the bias 
against finding KBOs in the 2:1 resonance over those in the 3:2 resonance. In §2, we set 
forth general, model-independent considerations for calculating this bias. In §3, we describe 
in detail the results of a particular simulation of resonant migration. In this section, we 
present illustrative snapshots of the instantaneous spatial distributions of Twotinos and of 
Plutinos. In §4, we explore how our results change by varying the migration rate of Neptune. 
In §5, we discuss our theoretical results in the context of the observations. There, we begin 
to examine critically the belief that Plutinos intrinsically outnumber Twotinos. A summary 
of our main findings is provided in §6. Our computations may serve not only to de-bias 
extant observations and thereby constrain the true relative resonant populations, but also 
to guide future observational surveys. 



The ability of an observational survey to detect KBOs residing within a given resonance 
depends on the KBOs' (1) spatial distribution, and (2) size and albedo distributions. In this 
section, we offer comments regarding the former consideration. 



By definition, an object inhabits a j + outer Neptunian resonance if the resonant 
argument, 



2. 



GENERAL CONSIDERATIONS 



2.1. 



Mean-Motion Resonances 



= U + 1)A- jAiV -u, 



(1) 
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librates (undergoes a bounded oscillation about a particular angle). Here j is a positive 
integer, A and Co are the mean longitude and longitude of periastron of the object, respectively, 
and Ajy is the mean longitude of Neptune. A restricted range for Qj+i-.j implies that the 
resonant particle will most likely be found at particular longitudes with respect to Neptune. 
For example, if an object inhabits the j '• = 1 resonance such that $2=1 librates about 180° 
with a negligibly small libration amplitude, then such an object attains perihelion when 
Neptune is 180° away in longitude. The eccentricity of the resonant object may be so large 
that the orbits of Neptune and of the particle cross, but the particle avoids close encounters 
with Neptune by virtue of the boundedness of $ 2: i- 

It has been remarked that because <3>3 : 2 for Pluto and the Plutinos librates about a mean 
value of ($3:2) = 180°, these objects tend to be found at longitudes displaced ±90° from 
Neptune when they reach perihelion and are at their brightest (e.g., Jewitt, Luu, & Trujillo 
1998). This argument is not strictly correct; it neglects the Plutinos' often substantial libra- 
tion amplitudes, A$ 3:2 . Just as a librating pendulum is most likely found near the turning 
points of its trajectory, a Plutino's resonant argument is most likely found near ($3:2) + A$3 : 2 
or near ($3:2) — A$3 : 2, not ($3:2)- Figures la and lb portray two toy models for the spatial 
distributions of 3:2 resonant objects. They demonstrate that the spatial distribution of reso- 
nant particles is sensitive to the distribution of libration amplitudes, dN/dA$, and not just 
to the value of the libration center, (<&). For each panel, the instantaneous locations of 15000 
co-planar particles are calculated according to the following scheme: semi-major axes are 
randomly chosen from a uniform distribution between 39.0 and 39.8 AU, mean longitudes 
are randomly drawn from a uniform distribution between 0° and 360°, and eccentricities are 
randomly selected from a uniform distribution between 0.1 and 0.3. Resonant arguments of 
particles are taken to equal §3-2 = ±180° + A^^sinA, where the upper and lower signs 
are equally probable and A is randomly drawn from a uniform distribution between 0° and 
360°. For Figure la, A$ 3 . 2 is randomly selected from a uniform distribution between 100° 
and 120°, and for Figure lb, the underlying distribution for A$ 3:2 is given by the solid his- 
togram in Figure 6. The longitude of perihelion of each particle is calculated according to 
Co = 3A — 2Ajv — $3:2, where \ N is assigned its present-day value of 302°. 

The resultant plots illustrate two ways that Plutinos could be distributed, both of which 
are possible in principle. In Figure la, the objects cluster in 4 locations, respecting the 2x2 
turning points of the resonant argument — 2 turning points for each of the 2 libration centers, 
($3:2) = ±180°. Thus, it is not true a priori that 3:2 resonant objects must cluster at only two 
locations in the sky. By contrast, in Figure lb, there are enough small-amplitude librators 
(A$ 3:2 < 1 rad) that the concentration of objects does gently peak at longitudes ±90° away 
from Neptune, in abeyance with the usual expectation. 
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Fig. 1. — Toy models for the spatial distribution of Plutinos. The particles' resonant ar- 
guments equal $3 : 2 = ±180° + A$3 :2 sinA, where A is uniformly distributed between 0° 
and 360°. In panel (a), A$ 3:2 is uniformly distributed between 100° and 120°. In panel 
(b), A$3 : 2 reflects the distribution obtained through simulation la, as shown by the solid 
histogram of Figure 6. Where Plutinos cluster depends sensitively on the distribution of 
A$ 3:2 . The dashed circles delimit radii of 40, 50, and 60 AU. 
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Toy models such as these are useful for analyzing the results of numerical orbit integra- 
tions. In anticipation of such integrations, we present Figure 2, which displays a toy model 
for the distribution of Twotinos. The parameters of the model are described in the figure 
caption. 

2.2. Secular Resonances 

Secular resonances might also play a role in determining the spatial distribution of 
KBOs. Chief among these are Kozai-type resonances in which u, the argument of perihelion, 
librates about particular angles, usually ±90°, 0°, or 180°. Pluto inhabits a Kozai resonance 
established by the total secular potential of all 4 giant planets, such that its u> librates 
about 90° with an amplitude of 23° (for a review, see Malhotra & Williams 1997). Thus, 
Pluto attains perihelion and is brightest only when it sits above the invariable plane by its 
orbital inclination of i p ~ 16°. If enough Plutinos inhabit Pluto-like Kozai resonances, their 
detection would be influenced by selection effects that depend on the latitude of observation. 

How many Plutinos and Twotinos inhabit Kozai-type resonances? Nesvorny, Roig, & 
Ferraz-Mello (2000) find that few observed Plutinos, ~2-4 out of 33, exhibit libration of 
oj. They explore a scenario by which Pluto gravitationally scatters other Plutinos out of 
the Kozai resonance. They find the scenario to be viable, though whether it is required by 
the model of resonant capture and migration is unknown; the efficiency of capture into a 
Pluto-like Kozai resonance in Malhotra's (1995) model of planetary migration might already 
be small enough to explain the observations. Unlike the case for the 3:2 resonance, we are 
not aware of any study of the possibility of calibration within the 2:1 resonance. In §3.4 
below, we investigate by direct numerical simulation the capture probability into Kozai-type 
resonances for both Plutinos and Twotinos within the model of planetary migration. 

3. MIGRATION MODEL 

Here we describe our model for the radial migration of the four giant planets and the 
resonant capture of planetesimals. Model ingredients are supplied in §3.1, mean-motion res- 
onance capture efficiencies are computed in §3.2, mean-motion resonance libration statistics 
and retainment efficiencies are discussed in §3.3, and the statistics of secular resonance cap- 
ture are presented in §3.4. Those readers interested in the spatial distribution of resonant 
objects may skip to §3.5 without much loss of continuity. 
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Fig. 2. — Toy models for the spatial distribution of Twotinos. Each panel portrays particles 
having libration centers, (<&2:i), that are indicated by the panel title. For ($2:i) = 180°, par- 
ticles avoid Neptune's longitude and are distributed symmetrically about the Sun-Neptune 
line. For ($2:i) = 90° (270°), particles cluster 90° ahead (behind) of Neptune's longitude. 
The combination of these libration centers in the proportions found in simulation lb yields 
panel (d), which resembles Figure 15. Toy model parameters are uniformly distributed over 
the following ranges: 47.2 AU < a < 48.4 AU, 0.1 < e < 0.3, 0° < A < 360°, 0° < A < 360°, 
135° < A$ 2:1 < 150° if ($2:1) = 180°, and 15° < A$ 2: i < 60° otherwise. 
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3.1. Initial Conditions and Migration Prescription 



To effect the migration, we follow Malhotra (1995) and introduce a perturbative accel- 
eration on each planet of the form 



where and a/ are the initial and final semi-major axes of a given planet, respectively, 
G is the gravitational constant, t measures time, r is a time constant, and v is the unit 
vector pointing in the instantaneous direction of the planet's velocity. Equation (2) corrects 
a typographical sign error in equation (7) of Malhotra (1995). This prescription causes each 
planet's semi-major axis to evolve according to 



but does not directly induce long-term changes in the planet's eccentricity and inclination. 
We adopt values for (aj, a/) for each of the planets as follows (in AUs): Jupiter (5.00, 5.20), 
Saturn (8.78,9.58), Uranus (16.2,19.2), and Neptune (23.1,30.1). 

We work in a coordinate system that takes the reference plane to be the invariable plane 
of the solar system. The positions and velocities of each planet at t — are adapted from 
Cohen, Hubbard, and Oesterwinter (1973), with the positions multiplied by a^/ay and the 
velocities multiplied by y/ 'af/ai. We employ the symplectic integrator, SyMBA, developed 
by Duncan, Levison, & Lee (1998), which is based on the algorithm by Wisdom & Holman 
(1991). The integrator was kindly supplied to us by E. Thommes, M. Duncan, & H. Levison. 

For simulations la and lb that are described in §3, we take r = 10 7 yr. Shorter migration 
periods of r = 10 6 yr and r = 10 5 yr are considered in §4. 

In simulation la, we focus on the efficiency of capture into, and the resultant dynamics 
within, the 3:2 resonance. We introduce 400 massless test particles whose initial semi-major 
axes lie between 31.4 AU (= 1 AU greater than the initial location of the 3:2 resonance) and 
38.5 AU (= 1 AU short of the final location of the 3:2 resonance). All particles in this region 
have the potential to be captured into the sweeping 3:2 resonance. Their initial eccentricities 
and inclinations are randomly drawn from uniform distributions between and 0.05, and 
between 0° and 1.4° = 0.025 rad, respectively. Arguments of periastron (cu), longitudes of 
ascending nodes (Q), and mean anomalies (M) are uniformly sampled between and 2n. 
The duration of the integration spans t 1 ? = 6 x 10 7 yr = 6r. 





a(t) = df — (cif — Oj) exp (—t/r) , 



(3) 
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In simulation lb, we concentrate on the 2:1 resonance. The only essential difference 
between simulations la and lb is that for the latter, the 400 test particles are distributed 
initially between 37.7 AU (= 1AU greater than the initial location of the 2:1 resonance) 
and 46.8 AU (= 1AU less than the final location of the 2:1 resonance). Thus, all such 
particles are potentially captured into the 2:1 resonance. The duration of this simulation is 
tf = 8 x 10 7 yr = 8r. 

3.2. Capture Efficiencies 

Of the 400 test particles in simulation la, 92 are captured into the 3:2 resonance. By 
definition, $ 3:2 librates for these 92 objects but circulates for the remaining 308. This 
capture efficiency of fs-_2 ~ 23% reflects (1) the probability of capture into the isolated 3:2 
resonant potential just prior to resonance encounter (Henrard & Lemaitre 1983, Borderies 
& Goldreich 1984), (2) losses due to pre-emptive capture into "competing" resonances such 
as the 5:3 and 2:1 which lie exterior to the 3:2, and (3) losses due to violent scattering by 
close encounters with the planets. Figure 3 displays the final eccentricities and inclinations 
versus the semi-major axes of those test particles having a < 60 AU. 

From simulation lb, we estimate the capture efficiency of the 2:1 resonance to be f 2: i ~ 
212/400 = 53%. This value is more than twice as high as / 3:2 , reflecting both the lack 
of competition from other sweeping resonances which lie interior to the 2:1, and the lower 
probability of scattering by Neptune at these greater distances. Figure 4 displays the final 
(a,e,i) for those test particles having a < 60 AU at the close of simulation lb. Note that 
many objects remain uncaptured by any low-order resonance at semi-major axes 43 AU < 
a < 47 AU; these non-resonant bodies presumably represent members of the low-inclination 
Classical Kuiper Belt that is observed today (Levison & Stern 2001; Brown 2001). 

These capture efficiencies are recorded in Table 1. We emphasize that / represents only 
the efficiency of capture, as distinct from the efficiency of retainment of captured objects, 
g, over the age of the solar system. Whether a captured KBO remains in a given resonance 
over 4 x 10 9 yr is discussed in §3.3. 

Mean inclinations, (i), of Plutinos and Twotinos are plotted against mean eccentricities, 
(e), in Figure 5. The mean is taken over the last 1 x 10 7 yr in simulation la, and the last 
3 x 10 7 yr in lb, during which times the migration has effectively stopped. There is a tendency 
for the Plutinos to have their (i)'s and (e)'s inversely correlated. No apparent correlation 
exists between (i) and (e) for the Twotinos. The predicted inclinations seem too low to 
compare favorably with the observed inclinations; the inadequacy of the migration model in 
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Final Eccentricities in Simulation la 
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Fig. 3. — (a) Final eccentricities vs. semi-major axes for particles with a < 60 AU in simu- 
lation la for which r = 10 7 yr. (b) Final inclinations vs. semi-major axes. Out of 400 test 
particles potentially swept into the 3:2 resonance, 92 are actually captured. Of these 92, 
perhaps only 42 would remain bound to the 3:2 resonance over the age of the solar system. 
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Fig. 4. — (a) Final eccentricities vs. semi-major axes for particles with a < 60 AU in simu- 
lation lb for which r = 10 7 yr. (b) Final inclinations vs. semi-major axes. Out of 400 test 
particles potentially swept into the 2:1 resonance, 212 are actually captured. Of these 212, 
perhaps ~100 would remain bound to the 2:1 resonance over the age of the solar system. 
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explaining the inclination distribution of Plutinos has been noted by Brown (2001). We note 
further that 4 out of the 5 observed Twotino candidates have orbital inclinations between 
11.8° and 13.5° — values characteristically larger than what the migration model predicts for 
objects in the 2:1 resonance. 

3.3. Libration Statistics and Retainment Efficiencies 

For the Plutinos in simulation la, we find that ($3:2) = 7r, as expected from resonant 
perturbation theory for low eccentricity orbits. The distribution of libration amplitudes, 
A$3 : 2 = [max($ 3 . 2 ) — min(<3>3 : 2)]/2, is supplied in Figure 6. Most objects have substantial 
libration amplitudes > 1 rad. Levison & Stern (1995) and Morbidelli (1997) calculate that 
Plutinos having large A$ can escape the 3:2 resonance over the age of the solar system. 
Thus, many of the Plutinos that are captured in our simulation la would not likely survive 
if we were to extend our integration to t 1 ^ = 4 x 10 9 yr. We define the retainment efficiency, 
#3:2, to be the fraction of captured Plutinos that either have A$ < 110° or that exhibit 
libration of 00 about ±90°. These selection criteria are motivated by the stability study of 
Levison & Duncan (1995; see their Figure 7). Of the 92 captured Plutinos in simulation la, 
42 satisfy our criteria for long-term residency; the resultant value for g 3:2 = 46% is recorded 
in Table 1. 

For the Twotinos in simulation lb, ($2=1) groups about 3 values: ~7r/2, 7r, and ~37r/2. 
The splitting of libration centers at large e from a single center at 71 to two additional centers 
at ~7r/2 and ~37i~/2 was explored by Morbidelli, Thomas, & Moons (1995), Malhotra (1996), 
and references therein. Figure 7 plots A$ 2: i against ($2=1) • Objects that reside more deeply 
in the resonance librate about ($2:i) ~ 37r/2 and 7r/2. The former center is more heavily 
populated than the latter at the level of 93 to 85 objects. Though this difference is formally 
statistically insignificant, we nonetheless believe that it reflects a heretofore unnoticed and 
physically significant signature of the migration model. Further evidence supporting our 
contention is provided in §4, where we demonstrate that faster migration speeds dramatically 
enhance the asymmetry to levels of statistical significance. Unlike for the case of the 3:2 
resonance, long-term, systematic studies of the stability of objects within the 2:1 resonance 
have not been published. It might be thought that the requirements for stability within the 
2:1 resonance are less stringent than for the 3:2 resonance, because the u s and z/ 18 secular 
resonances do not overlap the 2:1 (Morbidelli et al. 1995). Simulations by Malhotra (2002, 
personal communication) indicate that the retainment efficiency is roughly #2:1 = 50%, and 
we will adopt this value in this paper. 

Many of the objects having ($2:i) = [ m ax($ 2: i) + mm (^2:i)]/2 ~ 7r also librate from 
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Fig. 5. — (a) Average i vs. average e of Plutinos during the final 1 x f0 7 yr of simulation la, 
after the planets effectively cease to migrate, (b) Average i vs. average e of Twotinos during 
the final 3 x 10 7 yr of simulation lb. While (i) and (e) for Plutinos tend to be inversely 
correlated, no such tendency exists for Twotinos. 
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Libration Amplitudes in 3:2 
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Fig. 6. — Distribution of libration amplitudes of captured Plutinos in simulation la (solid; 
r = 10 7 yr), Ha (dotted; r = 10 6 yr), and Ilia (dashed; r = 10 5 yr). Simulations Ha and Ilia 
are discussed in §4. Reducing r below 10 7 yr yields smaller libration amplitudes and thus 
greater efficiencies of retainment of objects within the 3:2 resonance. 
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Fig. 7. — Amplitude of libration vs. libration center for Twotinos. The majority of Twotinos 
librate about ($2:i) ~ ±vr/2 and have libration amplitudes that are smaller than those 
of Twotinos librating about (<&2:i) = ^- The lobe at (<&2:i) ~ 3n/2 contains 10% more 
particles than the lobe at (<&2:i) ~ 7r /2; this asymmetry becomes more pronounced as shorter 
migration timescales are considered; compare with Figure 18. 
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time to time about ~7r/2 and ~37r/2. An example of such a "three-timing" Twotino is 
displayed in Figure 8. A histogram of libration amplitudes of Twotinos is given in Figure 9. 

For completeness, (e), (i), and A$j +1 .j are plotted against one another in 

Figures 10 and 11. One correlation that emerges is that between (e) and (<&2:i) for (<&2:i) ~ 
7r/2, 37r/2. Increasing (e) increases the separation of the two principal libration centers away 
from (<&2:i) — 7T, an effect seen previously by Malhotra (1996). Another correlation appears 
between A$ 3 . 2 and (i). Large amplitude librators tend to have large (i). 

3.4. Capture into Secular Resonances 

Of the 92 captured Plutinos in simulation la, 9-17 evince libration of u over the last 
1 x 10 7 yr of the simulation. Our uncertainty arises because for 8 Plutinos, the duration 
of the simulation is too short to see either a complete cycle of libration or of circulation. 
The libration centers, (u), are distributed over ~90°, ~180°, and ~270° for the 9 confirmed 
calibrators. A sampling of the time evolution of u for three calibrators in the 3:2 resonance 
is provided in Figure 12. 

Only a subset of the 92 captured Plutinos are likely to remain bound to the 3:2 resonance 
over 4 x 10 9 yr. Of the 92, we estimate that 42 are potential long-term residents: they either 
have A$3 : 2 < 110° or their cu's librate with small amplitude about ±90° [compare with 
Levison & Stern (1995)]. Of these 42, 8-12 are calibrators, or 20-30%. Since the calibrators 
comprise a minority among surviving 3:2 resonant objects, and since (cu) can equal 180° 
in addition to ±90°, we conclude that secular libration of cu resulting from the migration 
model probably does not introduce strong latitudinal selection effects for the discovery of 
Plutinos. Of course, a zeroth-order assessment of latitudinal selection effects for resonant 
objects probably requires that we move beyond, or in the extreme case abandon, the model for 
resonant migration, since it apparently cannot reproduce very well the observed inclination 
distributions; see §3.2 and Brown (2001). 

Similar conclusions are obtained for Twotinos. Of the 212 objects captured into the 2:1 
resonance, only 12-15 evince libration of u. Figure 13 samples three of them. 

Our results address a question posed by Nesvorny, Roig, & Ferraz-Mello (2000; see 
their footnote 3). Could the observed paucity of Plutinos in the Kozai resonance today be a 
primordial relic of resonant capture and migration? Of the 42 captured Plutinos in simulation 
la that might survive in the 3:2 resonance over the age of the solar system, only 3-5 occupy 
a Kozai resonance for which (uj) = 90°. This fraction of ~4/42 compares well with the 
observed fraction of ~3/33 reported by Nesvorny et al. (2000). Thus, our answer to their 
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Ib: A Three -Timing Twotino 
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Fig. 8. — A Twotino that alternates between librating about ($2:i) ~ 7r/2, it, and 37r/2. 
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Ib: Libration Amplitudes in 2:1 
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Fig. 9. — Distribution of libration amplitudes of Twotinos in simulation lb. 



-19- 



Simulation la 3:2 
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Fig. 10. — (a) Average e vs. average <3> for Plutinos. (b) Amplitude of libration vs. average 
e for Plutinos. (c,d) Same as (a,b) but for Twotinos. As seen in (c), Twotinos with higher 
eccentricities librate about centers closer to (<&2:i) — 0° than those with lower eccentricities. 
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Simulation la 3:2 



Simulation la 3:2 



CD 
T3 



15 


t 






180 




12 


• 
• 






150 


• • • 


9 


• 
* 




60 


120 


• • • * 






CD 




• •• . 




\ 
I 




T3 










CM 


90 


•••••• 


6 


\ 

1 




n 
O 
< 


60 


• • 


3 


t 






30 


* 




$ 
t 
t 


(a) 






(b) 





• 










) 90 180 270 360 
<* 3:2 ) (deg) 

Simulation lb 2:1 



) 3 6 9 12 15 
(i) (deg) 

Simulation lb 2:1 




QO 
CD 
T3 



■e- 
< 



180 
150 

120 

90 

60 

30 




• 




• 




• 




• • • 
"j " * * 


• 


(d) 



90 180 270 360 
<* a i> (deg) 



3 6 9 12 15 
(i) (deg) 



Fig. 11. — (a) Average i vs. average $ for Plutinos. (b) Amplitude of libration vs. average i 
for Plutinos. (c,d) Same as (a,b) but for Twotinos. As seen in (b), highly inclined Plutinos 
tend to be weakly bound to the first-order, eccentricity resonance. 
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la: cj-librators in 3:2 
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Fig. 12. — Time evolution of uo for three Plutinos that also inhabit Kozai-type resonances. 
Since the center of libration, (u), does not take a unique value, and since fewer than 30% 
of Plutinos are found in the simulation to be calibrators, we conclude that Kozai-type reso- 
nances do not introduce strong latitudinal biases for finding Plutinos. 
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lb: w-librators in 2:1 
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Fig. 13. — Time evolution of u for three Twotinos in Kozai-type resonances. Fewer than 7% 
of the Twotinos in simulation lb are w-librators. 
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question is yes — scattering effects by Pluto need not be invoked to explain today's observed 
low fraction of Plutinos in the Kozai resonance. We leave unaddressed two issues: (1) the 
effects of collisions amongst KBOs in the primordial belt in populating the Kozai resonance, 
and (2) whether the fact that the most massive Plutino yet discovered also inhabits the 
Kozai resonance is coincidental. 



3.5. Spatial Distribution of Resonant Objects 

We synthesize instantaneous snapshots of the spatial distribution of resonant objects 
from our simulation data as follows. Essentially, the positions of R resonant particles sampled 
at T different times are taken to represent the positions of R x T particles sampled at 1 
time. Take the 2:1 resonance as an example. The inertial Cartesian coordinates of Neptune 
[(xn, UN, zn)] and of the R = 212 Twotinos [(xi, Zi), i — {1, . . . , 212}] at a time, t, in the 
simulation are rotated about the z-axis by an angle 

9(f) = arccos + V»W» . (4) 

Here Xj^ = 30.1 cos (302°) AU and = 30.1 sin (302°) AU are the approximate current 
coordinates of Neptune. This procedure shifts the positions of all bodies into the Neptune- 
centric frame. We repeat this operation for T different instants of the simulation to generate 
T different snapshots. These snapshots are then overlaid on one another to yield a single 
image, a representation of the present-day spatial distribution of R x T resonant particles. 

As noted in §3.3, only a subset of captured Plutinos will be retained by the 3:2 resonance 
over 4 x 10 9 yr. We employ the subset of R = x / 3 . 2 x 400 = 42 objects to construct the 
Plutino snapshot plot (Figure 14) and those figures quantifying the longitudinal variations 
of Plutino density (Figures 16a and 17). 

Our method is not strictly justifiable if Neptune does not execute a perfectly circular 
orbit that remains in the invariable plane. In practice, Neptune's eccentricity and inclination 
are so small that we do not consider this a serious violation. A more weighty concern is 
whether the true distributions of orbital elements — eccentricities, inclinations, and libration 
amplitudes — of RxT particles are well represented by only R particles. Since R is not small 
for either simulation, and since the distribution functions displayed in Figures 6, 7, 9, 10, 
and 11 do not betray poor coverage of phase space, we proceed with confidence. 



Figures 14 and 15 showcase the present-day snapshots of Plutinos and Twotinos, re- 
spectively. For the former figure, T = 1000 time slices of R = 42 particles are sampled 
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uniformly between t± — 5 x 10 7 yr and t 2 = 6 x 10 7 yr; for the latter figure, T = 190, 
R = 212, ti = 5.000 x 10 7 yr, and t 2 = 5.190 x 10 7 yr. More time slices could be sampled 
for the Twotinos, but we restrict T to keep the plots legible. These snapshots resemble 
closely the toy models presented in Figures lb and 2d; we conclude that (weak) correlations 
between orbital elements such as A$2:i and (e) — correlations that are missing from the toy 
models — do not significantly influence the organization of resonant populations. 

Plutinos in the migration model cluster ±90° away from Neptune. Twotinos cluster 
±75° away. Plutinos avoid longitudes near Neptune and longitudes that are 180° away from 
Neptune. Twotinos largely avoid longitudes near Neptune. Longitudinal variations in the 
density of resonant objects are quantified in Figure 16, where we plot the number of objects 
per degree in longitude, Nj + i.j, inside a heliocentric distance, r. In constructing Figure 16, 
we employ all available time slices (T = 1000 in simulation la and T = 3000 in simulation lb), 
and then normalize the curves for Nj +i: j so that the total population within each resonance, 
integrated over all longitudes, equals 10000 objects. 

The sweet spots on the sky for finding resonant objects are sweetest — that is, the con- 
trast between maximum object density and minimum object density is greatest — for small 
limiting r (e.g., r < 40 AU). In the large r limit, the maximum contrast in object density is 
~200% for Plutinos and ~40% for Twotinos. Note that the two sweet spots for Twotinos dif- 
fer in strength; the spot displaced by —75° from Neptune contains more objects than the spot 
displaced by +75°, reflecting a greater population of objects librating about ($2:i) ~ 37r/2 
rather than ~7r/2. See §3.3 and §4 for more discussion of this asymmetry. 

In Figure 17, we divide N^ :2 by N 2 -i to compute the longitudinal bias in finding Plutinos 
over Twotinos. Though the population of each resonance is normalized to the same number, 
many more Plutinos will be found than Twotinos at small r, simply because the 3:2 resonance 
is located closer than the 2:1. Most interestingly, however, there exists a special longitude 
interval between 210° and 240° (between —90° and —60° of Neptune's longitude) where 
approximately equal numbers of Plutinos and Twotinos are expected to be found. Within 
this longitude interval, the yield of Twotinos to Plutinos ranges from 0.4 to 1 as the limiting 
r increases from 40 AU to oo (assuming that the two resonances are equally well populated). 
The absolute object density (Nj +1 .j) for each resonance is also maximal over this longitude 
range, making this interval the sweetest of spots. A similar spot exists between 0° and 40° 
longitude (+60°-100° of Neptune's longitude); here N 2: i/N 3:2 varies from 0.3 to 0.9 as the 
limiting r increases from 40 AU to oo. 
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Distribution of Plutinos in the Kuiper Belt 




x (AU) 

Fig. 14. — Synthesized snapshot, viewed from the invariable pole, of the spatial distribution 
of low libration amplitude Plutinos in simulation la, for which r = 10 7 yr. Large black 
dots mark the positions of observed Plutino candidates catalogued by the Minor Planet 
Center as of July 4, 2002. The main features of the snapshot, including the "sweet spot" 
concentrations of Plutinos displaced ±90° from Neptune's longitude, and the relative dearth 
of Plutinos at longitudes 0° and 180° from Neptune's, can be reproduced by a simple toy 
model; compare with Figure lb. Dashed circles indicate heliocentric distances of 40, 50, and 
60 AU. Radial lines delineate where the Galactic plane, ±10° Galactic latitude, intersects 
the invariable plane; Kuiper Belt surveys avoid the Galactic plane because fields there tend 
to be too crowded with stars. 
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Fig. 15. — Synthesized snapshot, viewed from the invariable pole, of the spatial distribution 
of Twotinos subsequent to resonant capture in simulation lb, for which r = 10 7 yr. Large 
black dots mark the positions of observed Twotino candidates catalogued by the Minor 
Planet Center as of July 4, 2002. The main features of the snapshot, including the "sweet 
spot" concentrations of Twotinos displaced ±75° from Neptune's longitude, and the relative 
dearth of Twotinos at longitudes near those of Neptune, can be reproduced by a simple 
toy model; compare with Figure 2d. Dashed circles indicate heliocentric distances of 40, 
50, and 60 AU, and radial lines delineate where the Galactic plane, ±10° Galactic latitude, 
intersects the invariable plane. Strangely, though the model predicts roughly equal numbers 
of Twotinos ahead and behind of Neptune's longitude, Twotino candidates have only been 
detected in the former lobe. Contrast this snapshot with the one in Figure 19, where we 
consider a migration phase that is 10 x shorter. 
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Fig. 16. — Longitudinal distribution of (a) Plutinos selected for likely long-term stability, 
and (b) Twotinos. Data are synthesized from simulation I, for which the migration timescale 
t = 10 7 yr. The number of objects within each resonance, integrated over all longitudes, is 
normalized to 10000. Each curve represents those objects located at heliocentric distances, 
r, less than the value shown. Adjacent curves differ by a multiplicative factor of 1.0955 in 
limiting r. Longitudes for Neptune and for the anti-Neptune direction are indicated. The 
contrast in object density from longitude to longitude is greatest at small limiting r. Plutinos 
cluster at longitudes ±90° away from Neptune, while Twotinos are most likely found ±75° 
away. 
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Simulation Ia + b / Normalized to 10000 
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Fig. 17. — Longitudinal variations of the bias in finding Twotinos vs. Plutinos in simulation I, 
for which r = 10 7 yr. As in Figure 16, the number of objects in each resonance is normalized 
to 10000. Panel (a) plots the fraction of Plutinos to Twotinos found inside a heliocentric 
distance, r, while panel (b) plots the inverse ratio. Adjacent curves differ by a factor of 
1.0955 in limiting r. At longitudes between 210° and 240° (between —90° and —60° of 
Neptune's longitude), yields of Twotinos and Plutinos are expected to be nearly equal. A 
similar interval exists between 0° and 40° longitude (+60°-100° of Neptune). 
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4. VARYING THE MIGRATION SPEED 

Simulations Ha and lib are identical to la and lb, respectively, except that we set 
t = 10 6 yr and t l J = 5 x 10 6 yr. For simulations Ilia and Illb, r = 10 5 yr and t 1 ^ 1 = 5 x 10 5 yr. 

Table 1 summarizes the computed capture efficiencies, /, and estimated retainment 
efficiencies, g, for each simulation. The retainment efficiency equals the fraction of captured 
objects that are expected to remain in the resonance over 4 x 10 9 yr. For Twotinos, we 
assume that g = 0.5 (see §3.3). For Plutinos, g equals the fraction of captured objects that 
either have A$ < 110° or that librate about (u) = ±90°. 

For r = 10 6 yr, / 2:1 ~ 15%, more than three times lower than the corresponding value 
for t = 10 7 yr. By contrast, /3 :2 (r = 10 6 yr) m 78%, more than three times as high as 
fz-.2{T = 10 7 yr) because fewer objects are lost to capture by the 2:1 or to close encounters 
with Neptune. For r = 10 5 yr, / 2:1 ~ 0%, while / 3:2 ~ 30%. Our results are consistent with 
those of Ida et al. (2000) and Friedland (2001). 

The Plutinos' distribution of libration amplitudes shifts substantially towards smaller 
values as r is reduced. Figure 6 plots histograms of A$ 3:2 for all three simulations, la, Ha, 
and Ilia. Table 1 records how the retainment efficiency, g, more than doubles for Plutinos 
as r is reduced from 10 7 yr to 10 6 yr, a consequence of the smaller libration amplitudes that 
characterize faster migration rates. 

Remarkably, objects fill the 2:1 resonance asymmetrically: captured Twotinos prefer to 
librate about (<E> 2: i) ~ 270° rather than (<& 2: i) ~ 90°. Figure 18 plots A$ 2: i against (<& 2: i) for 
simulation lib. The preferential filling of one resonance lobe over another dramatically affects 
the spatial distribution of Twotinos, as illustrated by Figure 19 and as quantified in Figure 
20. The difference in populations is at the level of 330% and is statistically significant. 2 
Establishing the relative populations of Twotinos ahead and behind of Neptune would offer 
a powerful constraint on the migration history of that planet. 

Figure 21 is appropriate for r = 10 6 yr and is analogous to Figure 17. 



2 We have verified that for r = 3x f 6 yr, the sign of the asymmetry remains the same and its magnitude is 
statistically significant and intermediate between that of simulations lb and lib — among the 141/400 objects 
captured into the 2:1 resonance, 3 x as many Twotinos librate about (<J>2 : i) « 270° than about ($2:i) ~ 90°. 
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Fig. 18. — Amplitude of libration vs. libration center for Twotinos in simulation lib, for 
which r = 10 6 yr. The lobe at ($2=1) — 270° boasts greater membership than the lobe at 
(^2:i) = 90° by a factor of 3.3-to-l. Contrast this result with Figure 7, for which r = 10 7 yr. 
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Twotinos Captured Quickly 
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Fig. 19. — Snapshot of the spatial distribution of Twotinos in simulation lib, viewed from 
the invariable pole. Large black dots represent observed Twotino candidates catalogued by 
the Minor Planet Center as of July 4, 2002. The shorter migration timescale of r = 10 6 yr 
leads to a pronounced asymmetry of objects in space. Dashed circles delimit heliocentric 
radii of 40, 50, and 60 AU, and radial lines indicate the intersection of the Galactic plane 
with the invariable plane. 
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Fig. 20. — Longitudinal distribution of (a) Plutinos and (b) Twotinos in simulation II for 
which t = 10 6 yr. The curves are normalized so that each resonance contains 10000 objects; 
adjacent curves differ by a factor of 1.0955 in limiting r. For Twotinos, the sweet spot 
located —75° away from Neptune's longitude is substantially sweeter than the one at +75°. 
Compare with Figure 16 which portrays the results of simulation I for which r = 10 7 yr. 
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Fig. 21. — Longitudinal variations of the bias in finding Twotinos vs. Plutinos in simulation 
11, for which r = 10 6 yr. As in Figure 20, the number of objects in each resonance is 
normalized to 10000. By contrast with simulation I for which r = 10 7 yr, only the interval 
between 210° and 240° longitude (—90° to —60° of Neptune's longitude) exhibits minimal 
bias in favor of detecting Plutinos over Twotinos; there, the yield of 3:2 resonant objects 
over 2:1 resonant objects ranges from 1 to 2.5 as the limiting heliocentric distance, r, out to 
which objects are detected increases from 40 AU to oo. 
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5. 



DISCUSSION 



Under the hypothesis of resonant capture, the number of Plutinos having diameters 
greater than s, divided by the number of similarly sized Twotinos, is given today by 



Here i] (> s, a) equals the number of objects having diameters greater than s that occupied 
an annulus of heliocentric radius a and radial width pa 8 AU just prior to the era of Neptune's 
migration. From our findings in this paper, (fz^Qz^) I {.f2:i92:i) equals 0.40, 9.6, and oo for 
r(yr) = 10 7 , 10 6 , and 10 5 , respectively (see Table 1). 

We interpret the observation that the 2:1 resonance today contains at least one (candi- 
date) object having a large eccentricity to imply that the migration timescale, r, cannot be 
as low as 10 5 yr. More refined estimates of r can be made by measuring the relative number 
of Twotinos whose resonant arguments ($2:i) librate about ~37r/2 rather than about ~7r/2; 
i.e., the relative number of Twotinos observed to reach perihelion at longitudes behind of as 
opposed to in front of Neptune's longitude. 

Estimates for the intrinsic F cannot be established without first de-biassing F f, s . As 
Figures 17 and equation (5) attest, F obs depends not only on r and the ratio of primeval 
populations, rj (a pa 35 AXJ)/r] (a pa 42 AU), but also on (1) the longitude of observation, (2) 
the limiting magnitude of the observation, (3) the relative size distributions of Twotinos and 
Plutinos, (4) the relative albedo distributions, (5) the relative inclination distributions, and, 
because of consideration (5), (6) the latitude of observation. For purposes of discussion, 
let us assume that Plutinos and Twotinos follow the same size, albedo, and inclination 
distributions, so that we can focus on considerations (1) and (2) exclusively. 

Roughly speaking, most KBOs have been detected by surveys having limiting magni- 
tudes my ~ 24 (Millis et al. 2002; see also the survey statistics compiled in Chiang & Brown 
1999). For this limiting magnitude, all objects that are inside r pa 44 AU and that have sizes 
s > 200 km and albedos > 0.04 would be detected. These sizes and albedos are comparable 
to those estimated by the Minor Planet Center. Moreover, a glance at Figures 14 and 15 
reveals that indeed, all but one of the observed Plutino and Twotino candidates have been 
discovered at r < 44 AU. 

Thus, if we assume that r = 10 7 yr, and take the dotted curve for r < 44 AU in Figure 17 
as our guide, then we crudely estimate the bias in finding Plutinos over Twotinos, averaged 
over all longitudes for which Plutinos have been discovered (330°-90°, and 150°-270°), to 




hajr) g 3 -.2 (t) t/ (> s, a pa 35 AU) 

/2:l(T) 02:1 (t) Vo(> S, d « 42 AU) ' 



(5) 
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be N 3:2 /N 2: i ~ 2.2. Then the observed ratio, F obs ~ 8, should be de-biassed down to values 
closer to 8/2.2 ~ 3.6. We emphasize that this is a model-dependent estimate of F that 
assumes that r = 10 7 yr. Since (/3 : 25'3:2)/(/2:ifi , 3:2)|T=io 7 yr ~ 0.40, we would estimate the 
ratio of primeval populations to be i] Q (a « 35AU)/?7o(a ~ 42 AU) ~ 3.6/0.40 ~ 9. This 
would reflect a steep drop in mass density over a short distance in the ancient planetesimal 
disk. Nonetheless, it would be consonant with the idea that a "Kuiper Cliff" (Chiang & 
Brown f999; see also Allen, Bernstein, & Malhotra 200f; Jewitt, Luu, & Trujillo 1998; 
Gladman et al. 1998) delineates the edge of the classical Kuiper Belt at a 48 AU. 

If we assume instead that r = 10 6 yr and repeat the same analysis by using the dotted 
line for r < 44 AU in Figure 21, we estimate an intrinsic F ~ 8/3.0 ~ 2.7, and a correspond- 
ing ratio of primeval populations of r] (a « 35A\J)/r] (a « 42 AU) ~ 2.7/9.6 ~ 0.3. Thus, 
while our estimate that Plutinos intrinsically outnumber Twotinos by a factor of ~3-to-l 
seems robust to changes in r, a dramatic drop in mass density with distance in the primordial 
planetesimal disk is by no means an assured conclusion. 

The above estimates for F are plagued by other observational biases that are difficult 
to quantify. For example, the Minor Planet Center dataset probably contains more Plutinos 
relative to Twotinos than it should because (1) MPC orbit fitting algorithms for objects dis- 
covered a few AU inside Neptune's orbit and having short astrometric arcs favor Plutino-like 
orbits as trial solutions (B. Marsden 2002, personal communication), and (2) the astromet- 
ric recovery rate of Plutinos is probably higher than that of Twotinos because the former 
objects are, on average, closer and therefore easier to re-detect by virtue of their brightness 
and large apparent proper motion. Both these factors lead us to conclude that our above 
estimates for F should be considered upper limits. 

Superior estimates for F can be obtained by coupling our theoretical calculations to 
the results of surveys having well-documented discovery statistics and minimal bias in their 
algorithms for orbit fitting. The Deep Ecliptic Survey (DES; Millis et al. 2002), for exam- 
ple, promises to be one such survey; it employs the more objective method of Bernstein 
& Khushalani (2000) in fitting orbits with short arcs. We defer analysis of the resonant 
populations using their large and homogeneous dataset to future study. 

An early but intriguing comparison between theory and observation lies in the complete 
absence of observed Twotino candidates at longitudes behind that of Neptune (see Figures 
15 and 19). By contrast, our numerical experiments demonstrate that resonant capture and 
migration preferentially fill the 2:1 resonant lobe displaced behind of, rather than ahead of, 
Neptune's longitude. Moreover, this asymmetry is only enhanced by faster rates of migration 
(Figure 19). With only ~5 candidate Twotinos, the probability of finding all of them in the 
forward lobe and not the backward lobe is 1/32, under the prior that a Twotino is as likely 
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to be found in one lobe as the other. Whether the actual Kuiper Belt defies the predicted 
sign of the asymmetry — in which case the present theory of resonant capture and migration 
must be considered either incomplete or incorrect — only continuing surveys for KBOs will 
tell. 



6. SUMMARY 



We have analyzed quantitatively the predictions of the model of resonant capture and 
migration for the 2:1 (Twotino) and 3:2 (Plutino) populations of the Kuiper Belt. We 
summarize our main findings as follows: 



1. The instantaneous spatial distribution of resonant objects depends not only on the 
libration centers of their resonant arguments, (<&), but also on their distribution of 
libration amplitudes, A$. For example, if A$ > lrad for most Plutinos, the usual 
expectation that such objects are most readily found ±90° away from Neptune's lon- 
gitude would not be valid. The distribution of libration amplitudes within a given 
resonance is model-dependent. 

2. We have numerically evaluated the capture efficiencies, /, of the sweeping 2:1 and 3:2 
resonances for three values of the migration timescale, r(yr) = 10 7 , 10 6 , and 10 5 . The 
timescale is assumed to be one of exponential decay. We define the capture efficiency 
to equal the fraction of objects whose orbits are initially spread uniformly over the 
complete path of a given sweeping resonance and that are ultimately captured by it. 
This efficiency depends not only on the probability of capture into the isolated resonant 
potential of interest, but also on the probability of pre-emptive capture into other 
resonances that lie exterior to the given resonance, and on the probability of violent 
scattering by close encounters with the planets. The capture efficiency for Twotinos, 
/2:i, decays from 53% to 0% as r is reduced from 10 7 yr to 10 5 yr. The capture efficiency 
for Plutinos, f^, increases from 23% to 78% as r decreases from 10 7 yr to 10 6 yr; the 
shorter migration timescale breeds fewer close encounters and fewer objects are lost to 
the competing 2:1 resonance. For r = 10 5 yr, / 3:2 ~ 30%. 

3. At least 1 Twotino candidate having a large eccentricity has been observed. This 
observation, interpreted within the confines of our migration model, implies that r 
cannot be equal to or lower than 10 5 yr. This conclusion is subject to the caveat that 
the migration model does not fully account for the observations; e.g., the model fails to 
generate the large orbital inclinations that are observed throughout the Kuiper Belt. 
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4. We have simulated the instantaneous spatial distributions of Twotinos and Plutinos 
as predicted by the migration model. If r = 10 7 yr, Twotinos cluster at longitudes 
displaced ±75° away from Neptune's longitude, where the upper sign corresponds to 
those objects librating with low amplitudes, A$ 2:1 < 1 rad, about (<&2:i) ~ 7r /2, and 
the lower sign corresponds to objects librating with similarly low amplitudes about 
(<&2:i) ~ 37r/2. Plutinos cluster at longitudes displaced ±90° away from Neptune's 
longitude, and all librate about ($3:2) = 7r. Longitude-to-longitude variations in the sky 
densities of Plutinos and Twotinos persist even for surveys that are not volume-limited 
in their ability to detect resonant objects of a given size. These variations sharpen as 
the limiting distance out to which resonant objects can be detected decreases. 

5. Over the longitude interval 210°-240° (—90° to —60° of Neptune's longitude), the bias 
in finding Plutinos over Twotinos is minimized. If the population of one resonance is 
identical to the other in terms of number, sizes, albedos, and orbital inclinations, then 
the migration model for r = 10 7 yr predicts 0.4 to 1 times as many Twotinos to be 
found over this longitude interval than Plutinos, as the limiting distance out to which 
objects are detected increases from 40 AU to 00. If r = 10 6 yr, the bias over this 
special longitude range varies from 0.4 to 1.1. A similar interval exists between 0° and 
40° (+60°-100° of Neptune's longitude) if r = 10 7 yr, but does not exist if r = 10 6 yr 
(see next point). 

6. The 2:1 resonance fills asymmetrically in the migration model — more objects are cap- 
tured into libration about ($2:1) ~ 37r/2 than about ($2:i) ~ 7r /2. The difference 
betwen populations is ~10% if r = 10 7 yr, and increases to ~330% if r = 10 6 yr. 
We reserve an analytic explanation for our numerical discovery to future study. 3 The 
asymmetry in libration centers translates directly into an asymmetry in the instan- 
taneous spatial distribution of Twotinos — more Twotinos are expected to be found 
at longitudes behind that of Neptune than in front of it. A differential measurement 
of the Twotino density ahead and behind of Neptune would powerfully constrain the 
migration history of that planet. 

7. Measuring the relative populations of Plutinos to Twotinos is a model-dependent en- 
terprise. Under the assumption that Twotinos and Plutinos share the same sizes, 
albedos, and orbital inclinations, we employ the results in this paper to crudely de- 
bias the current tally of ~43 observed Plutinos to ~5 observed Twotinos to estimate an 



3 We have performed additional numerical and analytic studies that demonstrate that this asymmetry can 
be understood purely within the restricted, circular, 3-body problem where the perturber's orbit expands 
outward. 
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intrinsic population ratio of F — 2.7-3.6. The range in our estimate reflects an order- 
of-magnitude difference in the assumed r, from 10 7 yr to 10 6 yr. We suspect that our 
estimate is an upper limit because Plutinos probably enjoy a greater frequency of astro- 
metric recovery than Twotinos, and because orbit fitting algorithms employed by the 
Minor Planet Center, whose dataset we use, favor Plutino-like trajectories. Improved 
estimates for F can be obtained by coupling our calculation to KBO surveys having 
well-documented discovery statistics. While it tentatively appears that Plutinos might 
intrinsically outnumber Twotinos by a factor not exceeding ~3, this conclusion can 
but does not necessarily imply that the surface density of the primordial planetesimal 
disk dropped dramatically with distance in the vicinity of ~42 AU. 
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Table 1. Summary of Simulations 



Label 


T(yr) 


Resonance 


/(%) a 


<?(%) b 


(*) f» 7T (%) C 


<*) » 3tt/2 (%) c 


(*) « tt/2 (%) c 


Figures 


la 


10 7 


3:2 


23 


46 


100 








3,5,6,10-12,14,16,17 


lb 


10 7 


2:1 


53 


50 


16 


44 


40 


4,5,7-11,13,15-17 


Ha 


10 6 


3:2 


78 


92 


100 








6,20,21 


lib 


10 6 


2:1 


15 


50 


13 


67 


20 


18,19-21 


Ilia 


10 5 


3:2 


30 


97 


100 








6 


Illta 


10 5 


2:1 
















a Efficicncy of capture into a given resonance. 

b Efnciency of retainment of captured objects by a given resonance over the age of the solar system. For the 2:1 resonance, 
this is assumed to be 50% (Malhotra 2002, personal communication). For the 3:2 resonance, we take the retained population 
to comprise those captured objects that either have A4>3 ; 2 < 110° or that inhabit Kozai-type resonances for which (w) = ±90° 
(compare with Levison & Stern 1995). 

c Percentage of objects in the resonance that librate approximately about the mean value indicated. 



